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Рассматривается задача вывода произвольной макроэкономической системы на магистральный путь развития по траектории, 
которая с точки зрения некоторого квадратичного критерия качества приближает пропорции валового внутреннего продукта к 
оптимальным пропорциям. При этом оптимальными пропорциями валового внутреннего продукта считаются пропорции маги- 
стральной (эталонной) экономической системы. 


Впервые П. Самуэлсон предложил именовать 
магистралью всякую сбалансированную траекто- 
рию с максимальным темпом роста. В связи с этим 
возникла бытовая трактовка достижения экономи- 
ческих целей оптимальным образом. А именно: для 
того, чтобы оптимальным путем «добраться» из на- 
чального пункта Х й в конечный пункт Х к следует 
быстрее выбраться на магистраль (улучшенную до- 
рогу для скоростного движения), «двигаться» по 
ней возможно дольше и лишь в конце свернуть в 
пункт назначения. 

Естественно, что перевод экономической систе- 
мы на магистральные темпы развития так же, как и 
переход к заданной цели, требует внешних усилий в 
виде инвестиционных вложений, причем возможно 
не разовых, а распределенных во времени. Для 
сложных многосекторных экономических систем 
от инвестора требуется еще и знание, в каких про- 
порциях инвестиции должны быть распределены в 
подсистемах или отраслях. Ответы на такие вопро- 
сы можно получить, рассматривая управление эко- 
номическими системами с позиций, развитых в 
теории оптимального управления и вариационного 
исчисления, теории систем и системного анализа. 

Любое управление предполагает наличие 
объекта управления (управляемой системы), 
субъекта управления (управляющей системы) и 
внешней среды. Основное назначение управляю- 
щей системы - поддерживать установленный и по 
каким-либо свойствам признанный нормальным 
режим работы объекта управления, а также обеспе- 
чивать нормальное функционирование отдельных 
элементов объекта управления в условиях воздей- 
ствия внешней среды. Объект управления во взаи- 
модействии с управляющей системой образует зам- 
кнутую систему управления. Связь, с помощью ко- 


торой управляющая система воздействует на 
объект управления, называется обратной связью. В 
теории автоматического регулирования функцию 
управления называют программой управления. 

В экономических задачах программа управле- 
ния соответствует принятию решения на перспек- 
тиву. Синтез оптимальных управлений в случае от- 
клонения состояния системы от планового значе- 
ния дает оптимальное решение и при новом ее со- 
стоянии. Изложенная ниже методика синтеза оп- 
тимальных управлений позволяет выводить дина- 
мические макросистемы на магистральные темпы 
функционирования . 

Постановка задачи. Имеем две динамические 
макросистемы, одну из которых будим называть 
развивающейся, 

Х(і) = В-\Е-А-К)Х(!) + В-'Щі)Х( 0)=Х0, (1) 
а вторую эталонной системой следующего вида: 

Х т {і) = С я Х т {1), Х т (0) = Х т 0, (2) 

здесь Х(і) и Х т (() - уровень валового внутреннего 
продукта развивающейся и магистральной систе- 
мы; А - матрица прямых материальных затрат; В - 
матрица приростных фондоемкостей; К - матрица 
конечного потребления; Е - единичная матрица; 
ІІ(() - внешнее инвестиционное воздействие; С т - 
матрица замкнутой магистральной системы. 

Системы дифференциальных уравнений (1) и 
(2) записаны в нормальной форме Коши и предста- 
вляют собой балансовые модели макроэкономиче- 
ских систем В. В. Леонтьева [1]. Моделирование и 
анализ динамических свойств таких систем, осно- 
ванный на контроле за расположением собствен- 
ных чисел на комплексной плоскости, был прове- 
ден Е.Л. Торопцевым [2]. 
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Особенности развития макроэкономики Рос- 
сии отражаются в системе национальных счетов 
(СНС). Неотъемлемой частью СНС являются та- 
блицы «Затраты-выпуск» [3], которые содержат по- 
дробные характеристики производства и использо- 
вания товаров и услуг, а также доходов, формирую- 
щихся в процессе производства. Симметричная та- 
блица «Затраты-выпуск», более известная среди 
российских пользователей статической информа- 
ции как межотраслевой баланс производства и ра- 
спределения товаров и услуг, устанавливает произ- 
водственные связи типа «отрасль-отрасль». В та- 
блицах собрана информация по 22-м отраслям. 

К сожалению, экономические параметры, пре- 
доставляемые ежегодным статистическим сборни- 
ком [3], можно использовать только в статических 
моделях; данная работа посвящена изучению дина- 
мических моделей. В качестве примера будут рас- 
смотрены малоразмерные трехотраслевые модели 
макросистем, которые получены путем агрегирова- 
ния. В качестве статических параметров для этих 
систем (матрицы А и К) использовались данные 
симметричной таблицы «Затраты-выпуск». Пара- 
метры системы, отвечающие за динамические ха- 
рактеристики модели (матрица В), а также началь- 
ные условия (2® и Х т 0), подобраны на основе про- 
ведения множественных имитационных вычисле- 
ний таким образом, чтобы была возможность на- 
глядной демонстрации переходных процессов, воз- 
никающих в макросистемах с инвестиционным 
воздействием. Экономические параметры моделей 
представлены в таблице. 

На рис. 1 представлены результаты интегриро- 
вания систем (1) и (2). Магистральное развитие 
эталонной системы (2) отражено на графике тон- 
кими линиями. В этой системе на всем временном 
интервале соблюдаются оптимальные пропорции 
развития, что приводит к постоянному росту вало- 
вого внутреннего продукта (ВВП). Для развиваю- 
щейся системы (1) внешнее инвестиционное воз- 
действие Щі)= 0, следствием чего является не опти- 


мальный и не сбалансированный план формирова- 
ния валовых выпусков, приводящий к угасанию 
одной из отраслей с течением времени. 


Таблица. Экономические параметры динамической модели 
развивающейся системы 


Прямые 
затраты (А) 

Конечное 

потребление 

(К) 

Приростные 

фондоемкости 

(В) 

Матрица 
замкнутой 
эталонной 
системы ( О т ) 

О 

О 

СО 

0,03 

О 

о 

сгт 

0,21 

0 

0 

0,31 

О 

о" 

0,67 

-2,48 

1,18 

1,52 

0,12 

0,23 

0,01 

0 

0,18 

0 

0,5 

0,32 

0,54 

1,7 

-1,81 

1,3 

0,02 

0 

0,09 

0 

0 

0,16 

0,48 

0,6 

0,2 

1,22 

1,04 

-2,33 


Развивающаяся система не в состоянии само- 
стоятельно перераспределить пропорции ВВП опти- 
мальным образом, поэтому в модели предусмотрена 
внешняя инвестиционная составляющая Щі), пока 
неизвестная, но благодаря которой должен полу- 
читься результат совмещения ВВП развивающейся и 
эталонной системы. Из этого следует, что разность 

г(і)=хм-х(о ( 3 ) 

стремиться к нулю при 1=со, т. е. 7(со)=0. Для прак- 
тических целей необходимо синтезировать такое 
внешнее управление 11(1), которое бы за конечное 
время 4 приводило разность 7(4) к нулю. При этом, 
начиная со времени 4, уровень ВВП развивающей- 
ся системы должен совпадать с уровнем эталонной 
системы. 

Решение задачи. Продифференцируем ур. (3) по 
1 и получим: 

У(і) = Х т (і)-Х(і). (4) 

Принимая во внимание (1) и (2), перепишем (4) 
в следующем виде: 

І(1) = С т Х т (1)-В-'(Е-А-К)Х(1)-В-'Ѵ(1), (5) 
учитывая, что Х т (і)=У(і)+Х(і), получим: 

У(1) = С т У(1) + 

+(С т -В-\Е- А- К))Х(1) - В-Щ). (6) 



Рис. 1 . Динамика ВВП магистральной и развивающейся системы без инвестиционного воздействия 
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В данном случае Х(1) совместно с Щі) являются 
управляющими функциями, с помощью которых си- 
стему (6) необходимо перевести в нулевое состояние. 
Определение 11(1) должно проводится в два этапа. На 
первом этапе необходимо приравнять (1(1) к нулю и, 
используя только управление от вектор-функции 
Х(1), привести систему (6) в устойчивое состояние. 
На втором этапе, зная траектории управления, мож- 
но определить 11(1), что будет показано ниже. 

Возможность перевода системы в нулевое со- 
стояние называется управляемостью или достижи- 
мостью. Для системы (6) матрица управляемости 
имеет вид [4]: 

Г с [а,В]=[р ар а 2 р ... а" ~'В\, (7) 

где а=С т , Р=С-В '(Е-А-К). 

Модель (6) полностью управляема тогда и толь- 
ко тогда, когда Г с имеет максимальный строчный 
ранг. Если система не полностью управляема 
(ранг(Г с )=А;<й), ее можно разложить на упра- 
вляемую и полностью неуправляемую подсистемы. 
В этом случае^ существует преобразование подобия 
Т, такое что У=Т 1 У, тогда: 

А = Т 'СТ- В = Т~\С т - В '(Е - А- К)), (8) 



’ А 

д. 


в ~ 

А = 

С 

12 

; в = 

с 


0 

Д,с_ 


0 


где А с имеет размерность к, и система (А с ,В е ) пол- 
ностью управляема. 

Управляемое подпространство модели (6) пред- 
ставляет собой любую возможную линейную комби- 
нацию состояний У с . Устойчивость этого подпро- 
странства определяется расположением собственных 
значений А с . Подпространство модели, не поддающе- 
еся управлению, состоит из всех состояний, получае- 
мых с помощью любой возможной линейной комби- 
нации состояний У К . Устойчивость этого подпро- 
странства определяется расположением собственных 
значений А к . Модель (6) стабилизируема, если ее под- 
пространство неуправляемых состояний устойчиво. 

Для упрощения дальнейших выкладок предпо- 
ложим, что модель (6) полностью управляема, тогда 
имеется возможность определения такого линейно- 
го квадратичного регулятора X, который удерживал 
бы выходы системы вблизи нулевого положения. 

Для замыкания системы введем следующее ли- 
нейное преобразование: 

Х(1)=2У(1). (10) 

Предположив, что 2 такой регулятор, который 
учитывает внешнее воздействие на систему (6) со 
стороны Х(і) и 11(1), получим замкнутую систему: 

У(1) = (С т + (С т -В\Е-А-К))2)У(1). (11) 

Для данной модели с желаемым набором соб- 
ственных чисел замкнутого контура можно вычи- 
слить набор постоянных коэффициентов усиле- 
ния, такой, что обратная связь по состоянию через 


постоянные коэффициенты усиления дает задан- 
ное расположение собственных чисел замкнутого 
контура на комплексной плоскости. 

С практической точки зрения реализация такого 
регулятора потребовала бы измерения величины 
каждой переменной состояния. Однако из-за эко- 
номических ограничений не всегда имеется воз- 
можность провести фактические измерения всех 
переменных состояния системы, в частности, ВВП 
отдельных отраслей. Это требует альтернативы фак- 
тическим измерениям и введения так называемых 
наблюдателей [4], иногда называемых также вирту- 
альными датчиками, фильтрами или расчетными 
данными. Цель наблюдателя состоит в том, чтобы 
вывести величину неизмеренной переменной со- 
стояния через другие переменные состояния, свя- 
занные с ней, и которые можно измерить. 

Наблюдатели имеют много общего с системами 
управления: 

• они являются динамическими системами; 

• с ними можно работать как в частотной, так и 
во временной области; 

• они могут быть проанализированы, синтезиро- 
ваны и спроектированы; 

• они имеют свойства такие, как устойчивость в 
переходных процессах и чувствительность; 

• эти свойства определяются полюсами и нулями 
их чувствительностей. 

В случае отсутствия полной информации о 
структуре ВВП эталонной и развивающейся систе- 
мы линейный наблюдатель для модели (11) будит 
следующим: 

т = (О я + (О я - В-\Е - А —К))2)У(і) + 

+У(СГ(0-СУ(0), (12) 

где матрица / называется усилением наблюдателя, 
С - матрица измерений, а У(1) - восстановленным 
состоянием. Заметим, что в частном случае, когда 
/=0, наблюдатель превращается в модель замкну- 
той системы (11). Именно этот случай нас будет 
интересовать в дальнейшем. 

Один из подходов к синтезу регуляторов с напе- 
ред заданными значениями собственных чисел 
замкнутого контура состоит в том, чтобы поме- 
стить динамические свойства системы в так назы- 
ваемый стоимостной функционал, который затем 
минимизируется количественно. Этот подход на- 
зывается оптимальным управлением, так как опти- 
мизируется функционал качества. Оптимизация 
переносит техническую задачу с непосредственно- 
го проектирования регулятора на проектирование 
критерия качества, из которого регулятор получа- 
ется автоматически. 

Определим 2 таким образом, чтобы использо- 
вание его в цепи положительной обратной связи 
(10) минимизировало квадратичный функционал: 
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лх) = ^у т ду+х т кх)сіі, (13) 

О 

здесь 0 - неотрицательно определенная, а К — по- 
ложительно определенная диагональная матрица 
весовых коэффициентов. Весовые матрицы 0 и К 
определяют соотношение между качеством регули- 
рования (как быстро процесс сходится к нулю) и 
затратами на управление. Функционал (13) являет- 
ся стандартным вспомогательным критерием по 
которому минимизируются отклонения между ма- 
гистралью и фактической динамикой ВВП разви- 
вающейся системы. 

Матрица 2 определяется путем численного ре- 
шения уравнения Риккати [4, 5] : 

2С„, +С Т Ш 2- 2(С т -В-'(Е-А-К))К-'х 

х{С т -В\Е-А-К)) Т 2+<д = 0. (14) 

Вычислив 2 и проинтегрировав систему (11), 
запишем ее решение в виде: 



Рис. 2. Динамика У(1) 


У(0 = е 1 -^- Ф ^'7(0). (15) 

Все собственные числа матрицы строго отрица- 
тельны, следовательно, процесс сближения траек- 
торий развивающейся и эталонной системы явля- 
ется устойчивым и с течением времени У(1) стре- 
мится к нулю. Данный результат отражен на рис. 2. 

Определение функциональной зависимости 
У(Г), которая оптимальна с точки зрения минимума 
функционала (13), позволяет вычислить траекто- 
рии сближения развивающейся и эталонной систе- 
мы по формуле: 

х„(0=хм-У((). (16) 

Результат вычисления (16) показан на рис. 3. 
Верхний пучок траекторий приближается к эталон- 
ному пучку до полного совпадения и продолжает 
свое развитие уже в магистральном режиме назна- 
ченных траекторий. 

Дальнейшее изложение посвящено количе- 
ственной оценке инвестиционного воздействия 


10 


Хт(і) 0 

Хт)!)-] 

Хт(1:)2 

Хпеѵѵ(1) 0 

Хпеѵѵ( I) і 

Хпеѵѵ(1) 2 

. 0.2 



Рис. 3. Сближение траекторий развивающейся и эталонной системы 
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Ѵ(і) на оптимально развивающуюся систему (16), 
траектории которой со временем приближаются к 
магистральным траекториям. 

Предположим, что инвестиционное воздей- 
ствие на развивающуюся систему должно быть 
пропорционально Х„(і), тогда: 

ЩІ)=-МХ П М, (17) 

где М- матрица инвестиционного воздействия, за- 
мыкающая систему (1) отрицательной обратной 
связью. 

Синтез М с использованием интегральных ква- 
дратичных оценок затруднен ввиду неустойчивости 
процесса расширения развивающейся экономиче- 
ской системы (1). Для этой системы квадратичный 
интеграл качества вида (13) будет расходящимся, 
что следует из требования постоянного расшире- 
ния ВВП, преследующего уровень магистральной 
траектории. Поэтому определение М проводилось 
по методу наименьших квадратов, минимизирую- 
щих функционал /'следующего вида: 

ПМ) (18) 

I 

В результате численной минимизации Р(М) бы- 
ла получена матрица М, коэффициенты которой 
характеризуют необходимый уровень инвести- 
ционных вложений на единицу ВВП: 



' 8,56 

-3,85 

-4,67" 


м = 

-5,93 

-2,34 

8,49 

, погрешность 


ѵ -3,68 

-3,16 

7,75 , 



аппроксимации Р(М)= 0,000005. 

Численная минимизация функционала (18) про- 
водилась с использованием квазиньтоновского мето- 
да программно-математической среды МаіСасІ. Про- 
верка гипотезы о нормальности распределения остат- 
ков была подтверждена на жестком 10 %-ном уровне 
по ^-критерию и критерию Колмогорова-Смирнова. 


Наличие положительных и отрицательных ко- 
эффициентов в М свидетельствует о двунаправлен- 
ности потока финансовых средств как в систему, 
так и из нее. Сравнительно большая величина этих 
коэффициентов говорит о больших финансовых 
потоках средств, необходимых для поддержания 
магистрального развития системы. 

Коэффициенты матрицы М позволяют опреде- 
лить временной график инвестиционных процес- 
сов Щі), рис. 4. Отрицательные значения вектор- 
функции 11(1) следует рассматривать как внешний 
инвестиционный поток средств, необходимый 
данной системе для магистрального функциониро- 
вания. Положительные значения 11(1) характеризу- 
ют экспортный финансовый поток, направленный 
из системы. 

Для данной макросистемы, взятой в качестве 
примера, на первоначальном этапе необходимы 
вложения в отрасль (7 0 ; при этом система будет ра- 
ботать в автономном режиме (без инвестиций) 
приблизительно 1 год. Далее потребуются инвести- 
ции в отрасль а затем почти одновременно для 
( 7 0 и \] г . Причем эти инвестиции должны постоянно 
увеличиваться, чего требует выбранное маги- 
стральное направление развития и соответствую- 
щие ему пропорции ВВП. 

В реальных экономических системах, например, 
макросистеме РФ, возможность полной управляемо- 
сти в рамках балансовой модели связана со значения- 
ми коэффициентов матрицы В приростных фондо- 
емкостей. Если В вырождена, что имеет место в слу- 
чае наличия всего двух фондосоздающих отраслей 
(«Строительство» и «Машиностроение») из 22-х, то 
возможности по управлению экономической дина- 
микой ограничены. Эти ограничения связаны с ма- 
лым количеством составляющих движения, уча- 
ствующих в переходном процессе. Для таких систем 
синтез матрицы М с использованием функционала 
(18) будет затруднен, что связано с трудностью обра- 
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щения вырожденной матрицы В. Поэтому реальную 
частично-управляемую систему необходимо разбить 
на полностью управляемую и неуправляемую части. 
При этом закон инвестиционного управления и ко- 
эффициенты матрицы М необходимо синтезировать 
в полностью управляемой части системы. 

Возможность определения коэффициентов ма- 
трицы М в замкнутой модели развивающейся си- 
стемы (1) свидетельствует о необходимости измене- 
ний в экономических параметрах модели для дости- 
жения поставленных целей магистрального функ- 
ционирования, так как без этой матрицы и соответ- 
ственно без инвестиционного воздействия система 
развивается несбалансированно. Сбалансирован- 
ность в системе может быть достигнута путем соот- 
ветствующих изменений коэффициентов матриц 
прямых затрат и приростных фондоемкостей, а так- 
же параметров, отвечающих за конечное потребле- 
ние, при этом матрица замкнутой системы (1) дол- 
жна оставаться неизменной. В этой матрице содер- 
жится информация о темпах роста и пропорциях 
магистрального развития каждой отрасли макроси- 
стемы. Изменение экономических параметров ди- 
намических систем непременно скажется на уров- 
нях финансово-денежных потоков, циркулирую- 
щих внутри системы и за ее пределами, но эти изме- 
нения можно контролировать с помощью матема- 
тического аппарата балансовых уравнений. 


СПИСОК ЛИТЕРАТУРЫ 

1. Леонтьев В. Межотраслевая экономика. - М.: Экономика, 
1997. - 478 с. 

2. Торопцев Е.Л. Моделирование процессов экономической ди- 
намики макросистем. - СПб.: Изд-во СПбГУЭФ, 2001. - 235 с. 

3. Соколин В.Л. Система таблиц «Затраты выпуск» России за 
2000 год: Стат. сб. / Госкомстат России. - М., 2003. - 116 с. 


Прогнозирование и поддержание соответ- 
ствующих уровней экспортиртно-импортных и ин- 
вестиционных потоков макросистем, проведение 
вычислительных экспериментов и выработка прак- 
тических рекомендаций входит в компетенцию де- 
партамента макроэкономического прогнозирова- 
ния Министерства экономического развития и 
торговли РФ. Определение конкретных сумм инве- 
стиционных потоков связано с дискретизацией не- 
прерывных моделей динамических систем. Возни- 
кающие при этом колебания и отклонения от за- 
данных траекторий необходимо будет компенсиро- 
вать дополнительными воздействиями. 

Вывод 

Предложенный математический аппарат про- 
ектирования оптимальных инвестиционных воз- 
действий на макроэкономические системы позво- 
ляет создавать и поддерживать магистральные тем- 
пы развития этих систем, проводить анализ и опти- 
мизацию временных инвестиционных процессов в 
фазовом пространстве отраслей, предоставляя при 
этом соответствующую информацию для инвесто- 
ров и лиц, принимающих решения в области мак- 
роэкономического прогнозирования, управления 
инвестициями и оптимального развития макроэко- 
номики. 
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